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ABSTRACT 


A  much  more  complete  and  detailed  picture  of  the 
external  environmental  y-ray  radiation  field  has  been 
obtained  through  y-ray  transport  calculations  of  exposure 
differential  energy  spectra,  integral  es^osure  rate 
spectra,  and  angular  exposure  rate  distributions  due  to 
sources  distributed  on  or  in  the  soil  half  space. 

The  radiation  field  is  examined  not  only  for  the  natural 
emitters  *  K,  U,  and  ®®®Th  but  also  for  y-rays  whose 
energies  are  typical  of  weapons  test  fallout.  The  energy 
spectra  and  exposure  rate  angular  distributions  are  shown 
to  vary  with  detector  height  and  source  distribution  with 
resulting  important  implications  in  regard  to  detector 
calibration  and  prediction  of  ground  level  exposure  rates 
from  aerial  survey  data. 

Exposure  rate  results  as  a  function  of  detector  height 
in  air  above  the  soil  air  interface  are  tabulated  for 
various  source  energies  and  source  distributions.  Although 
the  calculated  results  are  for  specific  soil  and  air 
densities,  soil  moisture  and  composition,  the  data  can 
easily  be  adapted  to  other  soil  and  air  conditions. 
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I.  INTRODUCTION 


A  knowledge  of  the  radiation  field  in  air  due  to  y-ray 
emitters  distributed  in  the  soil  is  important  for  estimating 
exposure  rates  from  a  given  source  concentration,  evaluating 
hazards  to  the  population,  and  properly  interpreting  and 
relating  radiation  measurements  made  at  ground  level  and  at 
airplane  altitudes.  The  sources  of  this  radiation  are  the 
naturally  occurring  radioisotopes,  fallout  from  nuclear 
weapons  tests,  or  unanticipated  releases  into  the  environment 
from  a  nuclear  reactor  installation. 

The  calculation  of  the  exposure  rate  in  air  due  to 
these  sources  in  the  soil  has  usually  been  based  on  infinite 
medium  buildup  factors.  The  exposure  rate  from  a  distributed 
source  in  the  soil  was  usually  calculated  by  assuming  a  single 
medium  composition,  either  all  air  if  the  sources  were  distri¬ 
buted  close  to  the  interface  or  all  soil  (aluminum  was 
usually  substituted  for  soil)  if  the  sources  were  distributed 
throughout  the  soil  half  space.  The  easily  computed  direct 
beam  exposure  rate  for  each  point  source  element  in  the 
chosen  medium  was  multiplied  by  the  appropriate  buildup 
factor  and  an  integration  performed  numerically  (or  analyt¬ 
ically  by  substituting  a  functional  fit  to  the  buildup 
factor  data)  to  obtain  the  total  esqjosure  rate^^^.  This 
rough  solution  was  necessary  because  of  the  lack  of  analytical 
transport  methods  for  solving  the  one-dimensional  two-media 
y-ray  transport  problem.  The  infinite  medium  buildup  factor 
approach,  besides  failing  to  account  for  the  differences  in 
the  transport  properties  of  soil  and  air  at  low  energies, 
provided  no  information  on  the  energy  and  angle  distribution 
of  the  Y-rays  entering  the  detector;  information  important 
for  interpreting  field  measurements. 

In  this  report,  we  present  accurate  calculations  of 
Y-ray  exposure  rates  in  air  (better  than  ±5%  S.D. )  for  source 
distributions  characteristic  of  natural  and  fallout  Y“^3y 
emitters  in  the  soil  as  well  as  differential  energy  spectra. 
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integral  exposure  rate  spectra,  and  exposure  rate  angular 
distributions  lieretofore  not  available.  In  addition,  we  discuss 
sensitivity  of  the  data  to  source  energy,  source  depth 
distribution  in  the  soil,  and  detector  height  above  the 
interface.  The  calculations  were  carried  out  using  a 
combination  of  the  P-3  and  DP-1  polynomial  expansion  matrix^ 
equation  method  for  solving  the  two  media  transport  problem 
Coii5)arison  of  these  results  with  those  obtained  using  the 
infinite  medium  moments  method  based  on  buildup  factors 
indicated  that  the  latter  method  was  probably  fairly  accurate 
(±10%)  for  detector  positions  close  to  the  interface  but  may 
have  been  significantly  in  error  far  from  the  interface, 
especially  for  low  energy  sources  distributed  in  the  soil 

In  addition  to  the  specific  calculations  of  the  differ¬ 
ential  energy  spectra,  integral  exposure  rate  spectra,  and 
exposure  rate  angular  distributions,  the  effects  of  varying 
the  soil  composition,  moisture  content,  and  density  are 
considered  together  with  an  evaluation  of  the  calculational 
accuracy. 
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II.  ADAPTATION  OF  THE  POLYNOMIAL  EXPANSION  MATRIX 
EQUATION  METHOD  TO  THE  AIR-SOIL  INTERFACE  PROBLEM 


Calculational  details  for  solving  the  two  media  transport 
problem  by  the  P-3  and  DP-1  polynomial  expansion  matrix 
equation  method  appear  in  a  previous  report^  In  general, 

the  method  consists  of  separating  the  spatial  and  angular 
dependencies  of  the  angular  flux  in  a  truncated  series  of 
orthogonal  polynomials  and  making  use  of  the  orthogonality 
properties  to  reduce  the  Boltzmann  equation  to  a  set  of 
coupled  integro-dif ferential  equations  for  each  spatial 
component  of  the  flux.  These  equations  are  solved  by  dividing 
the  energy  range  into  a  number  of  groups  and  replacing  the 
integration  over  energy  by  a  summ.ation  over  these  groups. 

The  integro-dif ferential  equations  are  rewritten  as  a  set  of 
differential  matrix  equations^®^  These  equations  have 
relatively  simple  exponential  solutions.  The  solution  for 
a  given  energy  group  constitutes  the  source  term  for  the 
next  lower  group.  Thus,  starting  at  the  source  energy,  the 
differential  energy  and  angle  spectrum  can  be  constructed 
stepwise  down  to  any  desired  energy.  Total  exposure  rates 
and  angular  exposure  rates  are  obtained  by  weighting  the 
differential  spectra  by  the  appropriate  energy  absorption 
coefficient  and  integrating  over  the  energy  and  angular 
intervals  of  interest.  All  cross  section  data  used  for 
these  calculations  were  taken  from  Hubbell^*^ 

The  calculations  discussed  in  this  report  are  for  a 
given  set  of  soil  and  air  densities,  soil  moisture  content, 
and  soil  composition.  However,  as  this  section  points  out, 
our  data  can  be  applied  to  other  soil  and  air  conditions. 

A.  Exposure  Rate  Dependence  on  Soil  Density 

The  calculated  exposure  rates  are  for  a  soil  in  situ 
density  of  1.6  gm./cm^  .  Actual  in  situ  densities  of  soils 
can  range  from  less  than  1,0  gm/cm^  to  over  2.0  gm/cm^ , 
although  a  typical  range^^^  would  be  1.1  to  1.8  gm/cm^ . 
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since  in  our  calculations  we  are  dealing  with  a  one¬ 
dimensional  medium  in  the  sense  that  only  the  direction 
perpendicular  to  the  soil-air  interface  enters  into  the  y- 
ray  transport  equation  and  then  only  in  terms  of  the  y- ray 
mean  free  path  (mfp) ,  the  effect  of  changing  the  soil 
density  is  equivalent  to  a  changing  scale  factor  in  this, 
direction.  This  is  because  the  y-ray  mean  free  path  is 
equal  to  the  inverse  of  the  total  attenuation  coefficient, 
a  quantity  directly  proportional  to  the  density  of  the 
medium.  Therefore,  the  soil  density  can  affect  calculations 
with  certain  source  distributions. 

For  a  source  distributed  exponentially  with  depth 
(S  =  where  S  is  the  activity  at  depth  z.  So  is  the 

surface  activity,  and  ot  is  the  reciprocal  of  the  relaxation 
length),  changing  the  density,  p,  is  equivalent  to  changing 
the  effective  source  distribution  such  that  a  =1.6  a/p  . 

(In  each  case  the  total  source  activity  in  a  1  cm  column 
=  1  y/sec) .  Thus,  increasing  the  soil  density  effectively 
buries  the  source  more  deeply,  decreasing  the  flux  and  expo¬ 
sure  rate  at  the  detector  in  the  air  half  space.  Therefore, 
in  order  to  use  the  data  given  at  certain  relaxation  lengths 
in  this  report  for  a  density  other  than  1.6  gm/cm  ,  one 
would  have  to  apply  the  above  transform,  a'  =  1.6  a/p  ,  to 
determine  which  relaxation  length  to  use.  Since  ct  “ 
corresponds  to  an  infinite  plane  source  it  follows  that  the 
soil  density  has  no  effect  for  this  case. 

When  our  calculations  are  for  a  uniformly  distributed 
source  of  one  y— ray  emitted  per  cm  at  1.6  gm/cm  ,  the 
calculated  exposures  and  fluxes  are  valid  for  a  uniformly 
distributed  source  of  intensity  1.6/p  gammas  emitted/cm  — 
sec  where  p  is  the  actual  Jji  situ  soil  density.  This  is  so 
because  changing  the  density  by  some  factor  is  equivalent  to 
changing  the  source  intensity  in  the  mean  free  path  interval 
(dt,  which  is  t  mfp  from  the  interface)  by  this  same  factor. 
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B.  Exposure  Rate  Dependence  on  Air  Density 


The  air  density  used  in  our  calculations  was  1,204  rag/cm^ 
corresponding  to  a  temperature  of  20°C.  Since  y-ray  dis¬ 
tances  must  be  measured  in  mean  free  paths,  the  fluxes  and 
exposures  for  other  air  densities  can  be  obtained  through 
adjustments  to  our  data  as  in  the  case  of  differing  soil 
densities.  Thus,  our  data  at  h  meters  represents  the  flux 
or  exposure  rate  at  a  distance  (1.204  x  10~^/p^)h  where 
is  the  new  air  density. 

C.  Exposure  Rate  Dependence  on  Soil  Composition  and  Moisture 
Content 


In  our  calculations,  the  assumed  soil  constituent  concen¬ 
trations  given  in  percentage  of  total  weight  were  SiOg  -  67.5%, 
AlgOa  -  13.5%,  FegOa  -  4.5%,  COg  -  4.5%,  and  Hg  0  -  10.0%. 

Changing  the  composition  of  the  soil  medium  affects  the 
fluxes  and  exposures  more  fundamentally  than  changes  in 
density  since  Compton  scattering,  which  is  the  dominant 
process  at  the  higher  energies,  is  proportional  to  the  average 
Z/A  of  the  soil  while  photoelectric  absorption,  which  dominates 
at  low  energies,  depends  on  the  atomic  number  Z.  Thus,  a 
soil  with  more  water  in  it  would  have  an  increased  Z/A  due 
to  the  lesser  number  of  electrons  per  hydrogen  atom  as 
compared  to  most  of  the  other  soil  constituents  and, 
therefore,  relatively  more  Conpton  scattering. 

High  Z  materials  in  the  soil  such  as  iron  could  also 
affect  the  exposure  rates  and  spectra.  We  investigated  this 
possible  dependence  on  soil  composition  by  comparing  calcu¬ 
lations  for  our  standard  mock  soil  which  has  an  average 
(Z/A)  of  ,503  with  similar  calculations  for  pure  aluminum 
having  a  (Z/A)  of  -.482,  a  range  of  (Z/A)  typical  of  most 
soils.  Relatively  minor  differences  were  found  in  the 
exposure  rates  calculated  for  these  two  materials  indicating 
that  the  radiation  field  is  not  extremely  sensitive  to  the 
exact  composition  of  the  soil.  For  a  ,662  MeV  plane  source 
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distributed  on  the  interface,  the  differences  for  both  the 
unscattered  and  total  exposure  rates  at  h  =  1  meter  were 
less  than  1%,  while  at  h  =  100  meters,  the  differences  were 
less  than  0.5%.  The  differences  were  slightly  larger  when 
the  source  was  distributed  in  the  soil.  For  a  uniformly 
distributed  1. 0  MeV  source,  the  unscattered  aluminum  exposure 
rate  was  found  to  be  about  4%  greater  than  the  corresponding 
soil  value  at  h  =  1  meter  and  about  3.5%  greater  at  h  =  100 
meters.  The  total  exposure  rates  were  3.4%  greater  for 
aluminum  than  for  soil  at  both  heights. 

A  further  comparison  was  made  by  calculating  the  exposure 
rates  for  soils  whose  ratios  of  constituents  other  than  water 
were  the  same  as  our  mock  soil,  but  whose  moisture  contents 
were  0%  and  25%  as  compared  to  the  10%  of  our  standard  soil. 

In  all  cases  the  density  was  kept  constant  at  1.6  gm/cm  and 
the  source  was  1.0  MeV,  uniformly  distributed,  since  the 
aluminum  calculations  indicated  the  greatest  differences 
occur  for  a  distributed  source.  The  0%  and  10%  moisture 
soils  gave  almost  identical  results  while  the  25%  moisture 
soil  gave  ejcposure  rates  only  1.7%  lower  than  the  10% 
moisture  soil  at  both  1  meter  and  100  meters.  This  small 
difference  is  in  the  direction  expected  since  the  increased 
moisture  content  would  cause  slightly  more  Coitpton  scattering 
relative  to  photoelectric  absorption.  The  (Z/A)  for  the  0% 
moisture  content  soil  was  .497,  while  the  (Z/A)  for  the  25% 
case  was  .512. 

Thus,  the  calculations  are  relatively  insensitive  to 
minor  differences  in  soil  composition  and  moisture  content 
and  should  be  valid  for  a  wide  range  of  soils.  Corrections 
would,  of  course,  have  to  be  made  when  the  density  changes 
due  to  rainfall.  These  corrections  are  discussed  under 
Section  II,  A. 

D.  Error  Estimates 


Exposure  rates,  in  this  report,  are  given  in  units  of 
MeV/gm-sec  and  represent  the  energy  absorbed  in  a  gram  of 
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air,  per  unit  time,  i.e.  the  flux  weighted  by  the  photon 
energy  and  mass  energy  transport  coefficient  for  air.  The 
flux  was  calculated  down  to  a  low  energy  cutoff  of 
To  +  13.2  ^Tq  =  This  low  energy  cutoff  has  been 

found  to  include  over  99%  of  the  energy  except  for  the 
highest  energies  and  for  very  large  interface  to  detector 
distances.  For  these  latter  cases  a  correction  was  made  by 
interpolating  the  differential  energy  spectrum  down  to  lower 
energies.  For  all  calculations,  annihilation  radiation, 
bremsstrahlung,  and  coherent  scattering  were  neglected. 

The  polynomial  expansion  method  can  be  expected  to 

provide  exposure  rate  estimates  for  scattered  y-rays  of 

accuracy  better  than  In  many  cases  the  scattered 

component  is  smaller  than  the  unscattered  component.  The 

latter  was  calculated  exactly  and  any  error  in  it  is  due 

only  to  errors  in  the  cross  section  data.  Thus,  the  error 

in  the  total  exposure  rate  is  smaller  than  the  error  in  the 

scattered  component.  To  improve  the  calculational  procedure, 

the  values  of  the  differential  scattered  flux  and  the 

esqjosure  rate  at  the  source  energy  were  calculated  directly 

using  an  exact  expression.  Inasmuch  as  the  y-ray  cross 

section  data  are  felt  to  be  quite  accurate  (better  than  2%) 

for  the  source  energies  and  media  used  in  our  calculations, 

we  conclude  that  the  error  in  exposure  rate  values  as  well 

as  in  the  differential  energy  spectra  and  the  integral 

exposure  rate  spectra  is  always  less  than  ±5%.  The  angular 

distributions  are  not  as  accurate  and  the  error  here  may  be 

as  much  as  ±10%  or  more  for  the  scattered  component  when 

(  2  ^ 

the  detector  is  near  the  interface 
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III.  RADIATION  IN  AIR  DUE  TO  THE  UNIFORMLY  DISTRIBUTED 
NATURALLY  OCCURRING  SOURCES  IN  THE  SOIL 


Gamma  radiation  in  air  due  to  naturally  occurring 
isotopes  in  the  soil  is  usually  calculated  by  assuming 
uniform  distribution  in  the  soil.  The  only  natural  emitters 
of  consequence  are  *^^U,  and  ^^^Th.  The  latter  two 

isotopes  decay  by  means  of  a  long  chain  process  and  K 
decays  to  '*°A  with  the  emission  of  a  single  1.46  MeV  y-ray. 
For  the  ®®®U  series,  most  of  the  y-rays  are  emitted  by  RaC 
(®^‘*Bi)  and  RaB(®^‘‘Pb)  which  are  daughters  of  *®®Rn^  a 
gaseous  daughter  of  *^^Ra.  Because  it  is  a  g^s^  Rn  can 
emanate  from  soil  particles,  diffuse  through  the  soil  air  to 
the  surface  and  escape  into  the  atmosphere.  This  effect 
tends  to  reduce  the  y-ray  emission  in  the  soil  due  to  U 
by  as  much  as  50%,  thus  lessening  the  importance  of  the  U 
series  relative  to  the  other  two  emitters  •  The  Th 
series  contributes  a  large  number  of  y-rays  of  varying 
energies  from  several  different  daughter  isotopes.  Although 
one  of  these  emitters  is  a  gas,  thoron  (^®°Rn),  the  short 
half  life  of  this  isotope  (54.5  sec)  prevents  any  significant 
escape  and  subsequent  loss  of  soil  gamma  activity. 

Since  the  decay  schemes  of  ‘^‘^^Ra  and  *®^Th  (particularly 
the  relative  intensities  of  the  various  decay  levels)  are 
still  somewhat  uncertain,  we  give  in  Table  1  the  source 
spectra  and  intensities  used  for  the  calculations  in  this 
report.  These  energies  and  intensities  are  our  estimates 
based  on  the  best  available  data.  Some  of  the  weaker  energy 
lines  have  been  grouped  together  or  included  with  stronger 
lines  in  order  to  keep  the  total  number  of  source  energies 
manageable. 

Table  2  gives  the  exposure  rates  for  monoenergetic 
sources  ranging  from  0.25  to  2.75  MeV  vs.  detector^height 
so  that  if  necessary  the  reader  can  recalculate  ^  U  and 
^^^Th  e5q>osure  rates  for  other  source  spectra.  The  ^  U 
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2  3  2 

and  Th  exposure  rates  were  obtained  by  interpolating  from 
the  data  in  Table  2  for  the  source  energies  in  Table  1, 
^ultiplying  by  the  number  of  y-rays  emitted  per  and 

Th  disintegration  and  summing. 

Calculations  of  the  exposure  rates  for  uniformly 
distributed  sources  in  the  soil  are  discussed  and  evaluated 
together  with  the  differential  energy  spectra,  integral 
exposure  rate  spectra,  and  angular  exposure  rates. 

A.  Exposure  Rates 

The  total  exposure  rates  as  well  as  the  direct  beam  or 
unscattered  component  are  given  in  Table  2  as  functions  of 
source  energy  for  various  detector  heights.  Figure  1 
illustrates  the  variation  of  exposure  rates  with  height  for 
three  different  source  energies.  Note  that  the  total  expo¬ 
sure  rate  changes  very  little  in  the  first  10  meters  above 
the  interface  but  then  begins  to  drop  off  fairly  rapidly  with 
height  with  the  lower  energy  sources  falling  more  rapidly 
than  the  higher  energy  sources  due  in  large  part  to  the 
more  rapid  decrease  of  the  unscattered  component  for  low 
source  energies.  The  scattered  components  fall  off  less 
rapidly  than  the  total  exposure  rate.  Since  the  calculated 
exposure  rates  for  series,  and  the  ®^^Th  series  all 

showed  almost  exactly  the  same  variation  with  height  even 
though  their  source  spectra  varied  considerably  (Table  3),  a 
single  curve  is  given  in  Figure  1  for  the  variation  with 
height  of  the  natural  gamma  emitters. 

The  variation  with  height  of  the  natural  emitter 
exposure  rate,  given  in  Figure  1,  was  crudely  verified  by 
us  at  a  single  location  by  making  ionization  chamber 
measurements  at  h  1  meter  and  h  7  meters  above  a 
predominantly  (95%)  natural  y-radiation  field.  The 
measured  variation  in  exposure  rate  with  height  of  14% 
compares  reasonably  well  with  the  11%  reduction  predicted 
by  Figure  1,  and  is  within  the  experimental  error. 
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Our  previous  calculations  of  Th  and  K 

exposure  rates  based  on  the  infinite  medium  buildup  factors 
are  also  given  in  Table  3.  These  calculations  for  an 
infinite  medium  assumed  that  the  energy  absorbed  at  the 
earth-air  interface  per  gram  of  air  would  be  5  the  energy 
emitted  per  gram  of  soil^^^.  Since  for  low  energies  the 
absorption  in  soil  is  greater  than  in  air  this  tended  to 
overestimate  the  exposure  rate.  By  comparing  the  K 
results  (Table  3),  we  see  that  this  overestimate  was  only 
about  4%  for  1.46  MeV.  It  would  have  been  greater  for  a 
lower  energy  source.  The  much  larger  differences  in  the 
old  and  new  calculated  exposure  rates,  however,  are  due 

not  as  much  to  the  more  accurate  method  (which  accounts  for 
only  about  5-8%  of  the  difference)  as  to  the  different 
source  relative  intensities  used  in  the  3  2 

any  case,  the  small  changes  in  the  values  of  U,  Th 
and  '*°K  exposure  rates  per  unit  concentration  of  these 
elements  do  not  significantly  alter  any  of  the  results 
given  in  any  of  our  previous  reports  on  environmental 
r adiat ion^ i • s ) _ 

The  present  report,  of  course,  is  much  more  detailed 
than  our  previous  work,  since  it  provides  data  on  energy 
and  angle  distributions  and  detector  height  dependence 
rather  than  just  exposure  rates  at  one  meter  above  the 
interface. 

B.  Differential  Energy  Spectra 


The  relative  contributions  of  various  energy  photons 
to  the  total  scattered  energy  spectrum  is  determined  by 
examining  the  differential  energy  spectra.  The  differential 
energy  spectra  of  the  scattered  energy  flux  (flux  x  energy) 
for  three  source  energies,  .364  MeV,  1.0  MeV,  and  2.5  MeV 
are  shown  in  Figure  2  for  h  -  1.0  meter  and  h  =  100  meters. 
The  effect  of  the  increased  scattering  in  air  relative  to 
soil  at  the  lower  energies  is  illustrated  in  the  figure  by 
the  shift  to  lower  energy  and  buildup  of  the  Compton  peak 


at  around  50-100  keV  as  the  detector  height  is  increased. 

The  relatively  high  photoelectric  cross  section  in  soil 
compared  to  air  causes  it  to  act  as  a  sink  for  low  energy 
y— rays,  an  effect  which  is  enhanced  the  lower  the  source 
energy.  To  further  illustrate  the  magnitude  of  this 
effect,  we  have  also  shown  in  Figure  2  the  spectrum  for  a 
1  MeV  source  at  1  meter  in  an  infinite  air  medium  instead 
of  a  soil-air  medium.  For  this  source  energy  the  difference 
in  the  total  exposure  rate  is  only  about  5%  in  all  since 
about  50%  of  the  exposure  rate  is  due  to  unscattered  y-rays. 
At  lower  source  energies,  where  the  scattered  y-rays 
contribute  a  larger  portion  of  the  total  exposure  rate,  the 
effect  on  the  total  exposure  rate  of  using  infinite  air 
calculations  would  be  more  significant. 


Integral  Exposure  Rate  Spectra 


Another  way  to  examine  the  effect  of  different  regions 
of  the  energy  spectrum  on  the  total  exposure  rate  is  to 
examine  the  integral  exposure  rate  spectrum,  i.e.  the 
fraction  of  the  total  exposure  rate  due  to  y-rays  of  energy 
less  than  E.  This  approach  can  be  extended  easily  to  the 
analysis  of  the  composite  and  ^®*Th  spectra.  Figure  3 

illustrates  that  as  the  detector  height  increases  the 
fraction  of  the  total  exposure  rate  due  to  low  energy 
photons  increases.  The  effect  is  more  pronounced  the  lower 
the  initial  source  energy  since  the  percentage  of  unscattered 
y-rays  is  lower  and  the  fractional  change  in  energy  due  to 
a  collision  is  smaller.  For  the  2.5  MeV  source  only  2%  of 
the  exposure  rate  is  from  y-rays  of  less  than  100  keV  at 
h  =  1  meter  and  only  5%  at  h  =  100  meters,,  while  for  a  .364 
MeV  source,  the  corresponding  values  are  13%  and  30%.  The 
values  at  the  source  energy  in  Figure  3  indicate  the 
fraction  of  the  exposure  rate  due  to- scattered  y-rays.  This 
fraction,  of  course,  decreases  with  increasing  energy. 

Figure  4  shows  the  integral  spectra  for  '*°K,  series 

and  ’“^‘“Th  series.  The  U  and  Th  spectra  are  quite  similar. 
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The  curves  which  represent  the  results  at  100  meters 
illustrate  the  softening  of  the  Spectral  composition  of  the 
radiation  with  height.  The  spectrum  at  1  meter  is  quite 

a  bit  harder  than  the  and  *®®Th  spectra  but  as  the 

detector  height  is  increased  this  difference  becomes  smaller 
since  the  higher  energy  sources  in  the  U  and  Th 
series  dominate  resulting  in  a  lesser  overall  rate  of  soften¬ 
ing  as  compared  to  that  for  the  single  1.46  MeV  '*  K  source. 
Figure  4  indicates,  however,  that  the  fraction  of  the 
exposure  rate  due  to  photons  of  less  than  200  keV  can  become 
quite  significant  as  the  detector  height  is  increased, 
depending  on  the  source.  Since  many  detectors  used  in  the 
field  have  an  energy  dependence  different  from  that  of  air 
at  low  energies,  it  may  be  inappropriate  to  use  the  same 
calibration  factors  for  these  instruments  at  different 
heights  above  the  interface.  This  would  be  especially  true 
for  an  instrument  calibrated  in  terms  of  photon  number  per 
unit  time  such  as  a  scintillation  counter  since  a  large 
increase  in  the  exposure  rate  due  to  low  energy  photons 
would  result  in  an  even  larger  increase  in  the  actual  number 
of  photons  at  that  height.  (This  can  be  seen  by  dividing 
the  points  in  Figure  2  by  energy  to  obtain  the  differential 
number  flux  curves) . 

The  softening  of  the  natural  emitter  y-ray  energy 
spectrum  with  detector  height  has  been  experimentally 
verified  qualitatively  in  the  field  both  by  Gustafson  et  al^^^ 
and  by  ourselves  in  similar  field  experiments  carried  out  in 
1965  using  Nal(Tl)  detectors. 

D.  Differential  Angular  Exposure  Rates 

The  differential  angular  exposure  rates  at  the  detector 
for  three  different  source  energies  at  h  =  1  meter,  100 
meters,  and  300  meters  are  shown  in  Figure  5.  These  curves 
are  given  in  terms  of  exposure  rate  per  radian  due  to 
photons  traveling  in  the  direction  0,  normalized  to  a  total 
exposure  rate  of  1.0,  where  9  is  the  angle  relative  to  the 
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perpendicular,  i.e.  0°  is  on  the  perpendicular  to  the  inter¬ 
face  heading  away  from  the  soil  half  space.  We  see  from  the 
figure  that  these  distributions  peak  at  around  70*’  -  80“  for 
h  =  1  meter,  40°  -  60°  for  h  =  100  meters,  and  30°  -  40°  for 
h  =  300  meters.  The  height  of  the  peak  decreases  as  the 
source  energy  decreases  while  the  "skyshine"  (photons  travel¬ 
ing  toward  the  ground)  increases  (see  Table  4) .  As  the 
detector  height  increases,  the  "skyshine"  for  a  given  source 
energy  also  increases.  The  flattening  out  of  the  angular 
distribution  as  the  source  energy  decreases  becomes  more 
pronounced  as  the  detector  height  is  increased.  The  angular 
distributions  of  the  scattered  and  unscattered  y-rays  are 
quite  comparable  (see  Figure  6)  to  the  total  exposure  angular 
distributions  although  the  unscattered  distributions  are 
slightly  more  peaked. 

The  angular  exposure  rate  distributions  for  '*°K, 
and  ®^*Th  (Figure  7)  are  similar  to  those  for  :364,  1.0,  and 
2.5  MeV  sources  shown  in  Figure  5.  There  is  little  difference 
in  the  three  distributions  for  h  =  1  meter  and  only  slightly 
more  for  h  =  100  meters.  Skyshine  contributes  11%  of  the 
*°K  exposure  rate,  13%  of  the  exposure  rate,  and  12%  of 

the  ®^*Th  exposure  rate  at  h  =  1  meter  and  12%,  15%,  and  13% 
of  the  exposure  rate  at  1  =  100  meters. 

The  natural  exposure  rate  angular  distribution,  therefore, 
is  fairly  insensitive  to  the  relative  amounts  of  ^°K, 
and  *^*Th  in  the  soil  although  it  does  vary  with  detector 
height.  Thus,  a  detector  with  an  angular  response  must  be 
calibrated  properly  with  elevation  in  order  to  interpret 
readings  of  natural  gamma  exposure  rates  made  at  various 
heights  above  the  interface. 
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IV.  RADIATION  IN  AIR  DUE  TO  EXPONENTIALLY 
DISTRIBUTED  (FALLOUT)  SOURCES  IN  THE  SOIL 


Radioisotopes  in  ttie  soil  as  a  result  of  fallout  from 
weapons  tests  (or  accidental  escape  from  reactors)  usually 
remain  on  or  just  below  the  surface  of  the  ground.  For  very 
short  times  after  deposition^  the  sources  can  even  be  assumed 
to  be  distributed  as  an  infinitesimally  thin  plane  source 
directly  on  the  interface.  For  longer  lived  fallout  which 
has  been  in  the  biosphere  for  some  time,  we  have  found  that 
a  reasonable  approximation  is  to  assume  the  sources  to  be 
distributed  exponentially  with  depth  in  the  soil  according 
to  the  relation  S  =  So  where  S  is  the  activity  at  depth 

z,  So  is  the  surface  activity,  and  a  is  the  reciprocal  of  the 
relaxation  length  (a  “  is  equivalent  to  a  plane  source  on 
the  interface) .  For  fallout  deposited  in  the  U.  S.  during 
1962  -  1965  an  assumed  relaxation  length  of  3  cm  gave  a 
reasonable  fit  to  actual  depth  distribution  measurements 

Many  investigators  calculate  the  exposure  rate  from 
fallout  isotopes  as  if  the  source  were  a  plane  source  buried 
at  some  depth  z  beneath  the  surface  in  order  to  account  for 
the  effect  of  ground  "roughness".  For  a  widely  distributed 
source,  for  detector  heights  at  least  a  meter  above  the 
interface  and  for  a  moderately  flat  surface,  the  assumption 
of  an  average  exponential  depth  distribution  should  be  more 
realistic  and  should  also  account  for  ground  roughness  since 
an  exponential  source  is  equivalent  to  burying  plane  sources 
successively  deeper  in  the  ground  with  decreasing  intensities 
The  choice  of  the  best  relaxation  length  to  use  for  a  given 
source  depends  on  such  factors  as  the  type  of  surface,  the 
source  energy,  the  time  since  deposition,  and  the  type  of 
soil. 


Variations  in  the  source  energy,  relaxation  length,  and 
detector  height  affect  the  exposure  rate,  energy  spectra, 
and  angular  distributions  in  air.  It  is  important  to 
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understand  the  significance  of  these  variations  to  properly 
interpret  as  well  as  predict  the  results  of  fallout  measure¬ 
ments  made  either  at  ground  level  or  from  an  airplane  or 
helicopter.  Thus,  the  results  of  our  confutations  of 
radiation  due  to  exponentially  distributed  sources  are 
discussed  in  this  light. 

A.  Exposure  Rates 

Total  ejqjosure  rates  for  exponentially  distributed 
sources  of  source  strength  one  gamma  emitted  per  cm®  of 
interface  surface  are  listed  in  Table  5,  A.  Unscattered 
ejfosure  rates  are  tabulated  separately  in  Table  5,  B.  The 
data  in  Table  5  should  be  sufficient  to  allow  the  reader  to 
construct  exposure  rate  vs.  height  curves  for  any  source 
energy  and  relaxation  length.  The  dependence  of  the 
exposure  rate  on  detector  height  is  shown  in  Figures  8,  9, 
and  10  for  several  source  energies  and  depth  distributions, 
(a  “  corresponds  to  a  plane  source.  All  our  calculations 
of  the  scattered  component  for  a  plane  source  were  arrived 
at  using  a  =  10,000.  The  unscattered  component  was  calcu¬ 
lated  exactly. ) 

From  Figures  8  and  9  we  see  that  the  variation  with 
detector  height  is  relatively  insensitive  to  source  energy, 
especially  below  h  =  100  meters.  Figure  10  illustrates  the 
effect  of  the  depth  distribution  on  the  esfosure  rate  at 
various  detector  heights.  The  scattered  component  falls 
off  very  slowly  with  height  all  the  way  up  to  about 
h  =  30  meters,  but  as  the  depth  distribution  approaches 
a  plane  source  the  unscattered  component  causes  the  total 
exposure  rate  to  begin  to  drop  off  more  quickly  with 
height.  The  exposure  rates  for  various  depth  distributions 
all  tend  to  converge  at  higher  altitudes,  i.e.  the  effect 
of  the  source  depth  distribution  is  reduced.  This  is 
qualitatively  what  we  would  expect  when  using  an 
exponential  source  distribution  model  to  represent  ground 
roughness,  since  ground  roughness  effects  decrease  as  the 
detector  height  is  increased. 
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From  Figure  11  we  see  that  the  ratio  of  the  exposure 
rates  at  h  =  1  meter  even  for  two  sources  quite  far  apart  in 
energy  does  not  depend  on  the  depth  distributiono  This  is 
important  since  often  one  can  make  the  assumption  that  two 
isotopes  are  distributed  similarly  with  deptho  If  one  then 
has  an  experimental  measure  of  the  ratio  of  their  activities 
at  any  depth  one  can^  using  the  results  of  this  report^ 
estimate  the  ratios  of  their  exposure  rates.  If  an 
independent  measure  of  the  total  exposure  rate  can  then  be 
made,  a  fairly  complete  picture  of  the  radiation  field  can 
be  deduced. 

A  more  detailed  examination  of  the  variation  in  exposure 
rate  with  detector  height  is  shown  in  Figure  12  where  the 
ratios  of  the  exposure  rates  at  10,  100,  and  300  meters  to 
that  at  1  meter  are  plotted  for  different  source  distributions. 
Here  we  see  that  the  10  meter/1  meter  and  100  meter/1  meter 
ratios  do  not  vary  too  rapidly  with  source  energy.  Sinc?e 
over  a  wide  range  of  depth  distributions  this  ratio  changes 
by  only  a  factor  of  about  2,  one  can  make  reasonable 
estimates  of  the  exposure  rate  at  ground  level  using 
measurements  made  from  an  airplane  or  helicopter  even  if 
the  exact  source  spectrum  of  the  radiation  is  unknown.  This 
fact  could  be  of  importance  for  some  types  of  emergency 
radiological  surveying  procedures . 

Figure  13  indicates  the  ratio  of  the  exposure  rate  due 
to  scattered  y-rays  relative  to  the  total  exposure.  The 
shape  of  the  curves  for  h  --  1  meter  and  h  “  100  meters  is 
fairly  similar,  reflecting  corresponding  dependence  on 
source  energy.  The  percentage  of  the  scattered  component  to 
the  total,  as  determined  from  Figure  10,  however,  is  much 
more  dependent  on  the  source  depth  distribution  at  h  1 
meter  than  it  is  at  h  —  100  meters.  This  again  indicates 
the  increased  effect  of  the  extra  soil  cover  on  exposure 
rates  at  lower  detector  heights  relative  to  higher 
altitudes. 

The  fractions  of  the  total  exposure  rate  due  to 
"skyshine"  are  given  in  Table  6.  The  "skyshine"  is  quite 
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sensitive  to  both  source  energy  and  detector  height  and 
becomes  a  very  significant  portion  of  the  total  exposure 
rate  at  low  source  energies  and  higher  detector  positions. 

B.  Differential  Energy  Spectra 

The  effect  of  the  source  depth  distribution  on  the 
differential  energy  spectra  of  the  scattered  energy  flux  is 
shown  in  Figure  14  for  a  .662.  MeV  (^®^Cs)  source.  At 
h  =  1  meter  the  spectrum  for  the  distributed  source 
(a  =  .33)  is  softer  than  that  for  the  plane  source  vhile  at 
h  =  100  meters  there  is  little  difference  in  the  shapes  of 
the  two  spectra.  Thus,  it  would  not  be  possible  to  make 
inferences  about  the  depth  distribution  or  the  exact 
exposure  rate  at  ground  level  by  using  an  altitude  measure¬ 
ment  of  the  energy  spectra. 

C.  Integral  Exposure  Spectra 

The  softening  of  the  scattered  energy  flux  spectra  at 
h  =  1  meter  res-ults  in  a  corresponding  softening  in  the 
integral  exposure  rate  spectra  (Figures  15,  16,  17) 
especially  for  h  =  1  meter.  All  three  figures  show  the 
same  general  features.  The  diminished  influence  of  the 
source  depth  distribution  at  h  =  100  meters  is  evident. 

The  fraction  of  the  exposure  rate  due  to  y-rays  below  a 
given  energy,  E,  increases  and  the  total  scattered  conponent 
increases  as  the  source  energy  decreases,  comparable  to  the 
result  we  obtained  for  the  uniformly  distributed  sources. 
Again,  we  must  emphasize  the  need  to  properly  interpret 
energy  dependent  dosimeter  readings,  since  at  h  =  100  meters 
over  20%  of  the  exposure  rate  from  an  source  (.364 

MeV)  is  due  to  y-rays  of  energy  less  than  100  keV  as  opposed 
to  a  corresponding  6%  for  1.25  MeV  C°Co)  y-rays. 

D.  Angular  Distributions 

For  a  plane  source  or  even  a  source  distributed  slightly 
with  depth,  at  detector  heights  close  to  the  interface  the 
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largest  portion  of  the  ej^osure  rate  is  due  to  unscattered 
Y-rays  (Figure  13) .  The  direct  beam  or  unscattered  y-raY 
angular  distributions  at  such  source  depths  rise  to  a  very 
sharp  peak  in  the  angular  region  0  =  89  —  90  ^  and  this 

peak  dominates  the  scattered  component  and  results  in  a 
very  skewed  total  exposure  rate  angular  distribution 
(Figure  18a) .  The  shape  of  this  distribution  is  slightly 
dependent  on  source  energy  becoming  less  peaked  as  the 
source  energy  decreases.  The  sharpness  of  the  peak  is  only 
slightly  diminished  as  the  source  becomes  more  deeply 
distributed  in  the  ground  and  is  still  very  acute  even  for 
a  =  .33.  The  angular  distributions  flatten  out  considerably 
as  the  detector  height  is  increased  but  even  at  h  =  100 
meters  (Figures  18a  and  18b)  these  are  still  more  peaked 
than  the  corresponding  distributions  for  uniformly  distri¬ 
buted  sources  (Figures  5  and  6) .  The  peaks  of  the  curves 
for  a  =  .33  (Figure  18b)  are  slightly  reduced  and  shifted 
toward  the  vertical  (9  =  o'*)  compared  to  the  corresponding 
curves  for  the  plane  source.  This  fact  might  be  useful  in 
estimating  the  depth  distribution  of  a  given  isotope  from 
an  angular  distribution  measurement  taken  in  an  airplane  or 
helicopter. 

Since  the  response  of  almost  all  instruments  has  some 
angular  dependence,  our  observations  are  significant  for 
interpreting  field  measurements.  It  is  clear  from  the  data 
presented  here  that  careful  instrument  calibrations  are 
essential  in  order  to  properly  interpret  measurements  made 
at  different  detector  heights. 
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V.  CONCLUSIONS 


The  energy  spectra  and  angular  exposure  rate  distributions 
for  fallout  sources  differ  considerably  from  those  for  the 
natural  emitters  because  of  the  different  source  depth 
distributions.  The  magnitude  of  this  difference  depends  on 
the  exact  depth  distribution  of  the  fallout  emitters.  The 
computations  in  this  report  indicate  that  both  distributions 
also  vary  with  detector  height ^  emphasizing  the  necessity  of 
carefully  calibrating  detectors  used  for  measuring  gamma 
radiation  in  the  field. 

The  calculated  exposure  rates^  differential  energy 
spectra,  and  angular  exposure  rates  can  be  used  for  inter¬ 
preting  and  analyzing  the  results  of  land  and  aerial  surveys 
over  extensively  contaminated  areas  as  well  as  for  predicting 
the  results  of  such  surveys  for  a  known  level  of  contamination. 
In  certain  instances  experimental  data  taken  at  high  altitudes 
can  be  used  to  infer  information  about  the  radiation  field 
near  the  interface.  This  can  be  done  by  utilizing  the 
curves  given  in  this  report,  taking  into  account  the  source 
spectrum  and  depth  distribution  of  the  isotopes  contributing 
to  the  field. 
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TABLE  1 


SOURCE  SPECTRA  USED  FOR  ^®®U  and  ®®®Th  EXPOSURE 
RATE  CALCULATIONS 
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Photons  per 

Energy 

Photons  per 

(MeV) 
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(MeV) 
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EXPOSURE  RATES  -  UNIFORMLY  DISTRIBUTED  SOURCE 


r>fnnc>irMc>ioic>ic>JCNr-«.-fr-ir-i 
•  I  »  *  *  •  w  J-ww  wsiv!.  w 

CMinOJOIfOfM'OO'rOOOi-IOOO 
^r^OfOVDvDa)0>fDroOr>*ro^ 


O  I  »  *  J^w 

S  TfcomMU^voroinTfiOTfinatH 
Oro<TvOC'ir-4in(X)r^r^'Huif>0 

c^iro^r^C0»H*-iiHCMr^nforo'^ 


^  ^Ttdoiror-vOvDvOmr-iHfori 
f>ir'VDroinO'D*-ioo^ovoa>r'i 

^vDO^rHi-HOJCNronTfinminvo 


CMCMr-4i-li-li-lrHfH.iHi-HrHi-lr-<iH 

I  I  I  I  I  I  I  I  I  I  I  I  I  1 


fo  00  o  Tf  O' 

n^cMwcomfSOOvOroO'VOO'ro 

inoowiHrHrsinm'^tininvOvDr^ 


m  m  n 
in  (N  00 


r-f  in  CNj  r-.  fM 
VO  ro  O'  VO  O'  m 


OJCMrHrHrHHrHHt-IrHi-C 

I  I  I  I  .1  I  I  I  I  I  I 


r-4voooin»-Hr-<nTtO'OvovO(N(r«“ 

r*.OcMr>or^^ooovOtNO'n'X> 


ino'i-irHrM(Mro^Tj‘invovot^r«- 


C>*C>ltHrHtHrHrHfHr~li-trH 

I  I  I  I  I  I  I  I  t  I  I 


oO’'tiNOr^O'roin(NfoocMO'in 

ooronooo.r'iniHOr*-^t*-i^oo 


inO'i-ir-irM<Nfo<<tvninvDr«.r‘r' 


CNCNr-lfHrHiH»Hr-l 

I  I  I  I  I  I  I  I 


I  I  I  I  I 


O'OOr^vD^fCO^OOOOOO'fMO'VO 

O'vOmOOHOOvOfNi-IO'inmvOO 


inO'MfHojcNn^ininvOr^r-OO 


I  I  I  I  I  I  I  I 


I  I  I  I  I 


vDO'oooor-»HoocNi-Hin^r-incM 

r-4r'moOfHO'VDroCNJO'vOror'iH 


voo'rHr-HiNojm^inmvDr'-r'OO 


I  I  I  I  I  <  1  t  t  I  I  1  I  I 


nO'OOO'^i-iininoO'nrHGO 

fMOO^O'fHO'r-nOJO'O^COfH 


^i)cr^^,_^f^J^M^OTtlnvOvDr-^'00 


CNJ  fN  iH  t-t  iH  iH 

I  1  I  I  I  I 


r— I  »— I  t— I  tH  I— I  iH 

I  I  t  I  I  I  1 


^<NO*-io'in<Nvo\o*HO^roa' 

oiO'^O'rHO'r^mnjor^^oOfH 


vDO'-H.HCMCNro^tnvDvDr^r^CO 


^  CM 

invoovoin  invOvOOinocMi.i 
(NninvDf-ocMttr-ocNinvOr-. 

,H.-<*-tr-4CN(NfMrsjrM 


<u 

u 

to 

1 

< 

u 

B 

z± 

V 

rH 

u 

0) 

•  - 

a 

u 

rC 

\ 

u 

fH 

II 

O' 

A 

in 

-P 

VO 

O' 

c 

II 

u 

U 

0) 

m 

to 

a 

1 

B 

o 

O' 

u 

v 

d 

> 

0 

(U 

CO 

s 

23 


24 


EXPOSURE  RATES  IN  AIR  DUE  TO  ‘‘°Kj  ®^®U,  and  ®^®Th  IN  THE  SOIL  (ixrAr) 
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TABLE  4 

PERCENTAGE  OF  TOTAL  EXPOSURE  RATE  DUE  TO  SKYSHINE  FOR 
UNIFORMLY  DISTRIBUTED  SOURCES 


Source 

ti  =  1  meter 

h  =  100  meters 

.364  MeV 

21% 

24% 

1.0  MeV 

13% 

16% 

2.5  MeV 

8.  6% 

9.5% 

®®®U  series 

13% 

15% 

‘‘®K 

11% 

12% 

®®*Th  series 

12% 

13% 

EXPOSURE  RATES  -  EXPONENTIALLY  DISTRIBUTED  SOURCE* 
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TABLE  5  (Cont'd) 
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TABLE  5  (Cont'd) 
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TABLE  5  (Cont'd) 
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TABLE  6 


PERCENTAGE  OF  TOTAL  EXPOSURE  RATE  DUE  TO  SKYSHINE- 
EXPONENTIALLY  DISTRIBUTED  SOURCES 


Source 

h  =  1  meter 

h  =  100  meters 

r  .364  MeV 

17.  6% 

23. 1% 

a  =  . 33  <  .662  MeV 

13.  9% 

18.  6% 

\  1. 25  MeV 

10. 1% 

13.3% 

r  . 364  MeV 

11.  8% 

23.1% 

Plane*  J  ^  662  MeV 

9.  5% 

H 

00 

s 

Source  1  ^25  MeV 

8.  0% 

12. 9% 

corresponds  to  an  infinite  plane  source. 


EXPOSURE  RATE  X  100  (MeV/gm-sec) 
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Figure  1.  Total  and  scattered  exposure  rates  versus  height  for 
uniformly  distributed  sources  with  energies  of  .364,  1.0  and  2.5 
MeV  and  the  single  normalized  curve  of  exposure  rates  for 
natural  emitters. 
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Differential  angular  exposure  rate  per  radian  (normalized  to  a  total 
rate  of  1, 0)  at  1^  100^  and  300  meters  for  uniformly  distributed  sources 


“H 

-H  C  C 
O  3 


ure  6,  Differential  angular  exposure  rate  per  radian  for  unscattered  gamma-rays 
y  (normalized  to  unscattered  exposure  rate  of  1.0)  at  1  and  100  meters 
formly  distributed  sources. 


. ^0  . 40  w  to  w  PPo  140  no 

9  (DEGREES) 

gure  7.  Differential  angular  exposure  rate  per  radian  for  natural  emitters 
rmalized  to  a  total  exposure  rate  of  1)  at  1  and  100  meters. 


HEIGHT  (METERS) 

Figure  8.  Total  exposure  rate  for  an  exponentially  distributed 
source  (a  .33]  for  various  source  energies  as  a  function  of 
height. 
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EXPOSURE  RATE  (MeV/gm-sec) 
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Figure  11.  Total  exposure  rate  as  a  function  of  energy  at  1  meter 
for  various  depth  distributions. 
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SOURCE  ENERGY  (MeV) 

Figure  12.  Effect  of  source  depth  distribution  on  total  exposure  rate  at  10^  100,  and 
300  meters,  relative  to  the  total  exposure  rate  at  1  meter. 


ENERGY  (MeV) 


Figure  14.  Differential  energy  spectra  of  the  scattered 
energy  flux  at  1  and  100  meters  for  a  . 662  MeV  source 
("^"Cs). 
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Figure  16.  Integral  exposure  rate  spectra  (fraction  of  total  exposure  rate  due  to 
gamma-rays  of  energy  less  than  E)  for  a  .662  MeV  source  at  1  and  100  meters. 
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igure  18b.  Differential  angular  exposure  rate  distributions  for  plane  and  exponentially 
istributed  sources  at  100  meters  (normalized  to  a  total  exposure  rate  of  1. 0) 


